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Abstract 

We present calculations on energy- and time-resolved two-photon photoemission spectra of im- 
ages states in Cu(100) and Cu(lll) surfaces. The surface is modeled by a ID effective potential 
and the states are propagated within a real-space, real-time method. To obtain the energy resolved 
spectra we employ a geometrical approach based on a subdivision of space into two regions. We 
treat electronic inelastic effects by taking into account the scattering rates calculated within a GW 
scheme. To get further insight into the decaying mechanism we have also studied the effect of the 
variation of the classical Hartree potential during the excitation. This effect turns out to be small. 
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FIG. 1: Schematic picture of a 2PPE experiment for a typical metallic surface with a surface band 
gap close to the Fermi level. The energy huj a corresponds to the pump photon and fiujb to the 
probe photon in 2PPE. On the right side we show the tail of the image potential, the densities 
of the surface state and first three image potential states of Cu(100) calculated with the model 
potential of Ref. and their corresponding binding energies. By Ep, E\ ower and E v&c we denote 
the Fermi energy, the lower edge of the energy gap and the vacuum energy. 



INTRODUCTION 



The presence of a metal surface creates electron states that do not exist in the bulk metal. 
For example, an electron situated at a distance z from the metal surface experiments an 
attractive force, F(z) = —e 2 / (2z) 2 , equivalent to that produced by its image charge situated 
at distance z inside the metal. For large z, the potential generated by this surface induced 
charge approaches the classical image potential, V(z) = -e 2 /Az^. When the metal has a 
surface band gap near the vacuum level, the electrons below the vacuum level are trapped 
between the well of the image potential and the surface barrier. These quantized states 
correspond to image potential states, and are spatially localized in the region in front of the 
surface (see Fig. [IJ. Metal image states typically have a large spatial extension, but barely 
penetrate the bulk metal. They form a Rydberg-like series with energies, E n , approximately 
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given by | 

En = -T^y, [eV] , (1) 

where n is a positive integer, and < a < 0.5 represents the quantum defect (which depends 
on the position and width of the gap). The energy of an electron trapped in an image state 
and moving parallel to the surface with momentum kii is therefore: 

E = —t + E n (2) 

2m e 

In contrast to bulk states, the small overlap between surface states and bulk states reduces 
considerably the inelastic electron scattering, leading to a long lifetime, t w of the image 
states, that scales asymptotically with the quantum number n as r n oc n 3 [2|. Due to this 
fact, these states are very interesting for the study of electron correlations. Furthermore, 
image states play an important role in the laser induced chemical control of reactivity at 
metal surfaces. It is therefore not surprising the considerable number of studies focusing on 
this subject (see Ref. [3] and references within). 

Image potential states of different surfaces have been observed experimentally^. A 
powerful technique to measure their energy and lifetime is two-photon photoemission spec- 
troscopy (2PPE)|5j. Recently, and with the advent of ultrafast laser technology it became 
possible to perform time- resolved 2PPE spectroscopy (TR2PPE), which allowed, with the 
technique of quantum beat spectroscopy jy, [D], the direct measurement of the lifetime of im- 
age states, even for states with large quantum number n. The 2PPE technique is depicted 
schematically in Fig. ^ One photon (pump), of energy hu a , excites an electron from an oc- 
cupied state below the Fermi energy (Ep) to an image-potential state of quantum number n; 
Then a second photon (probe), with energy tkui,, ejects the electron out of the surface, above 
the vacuum energy. This electron has a kinetic energy Ek = hu>b — E n , that is measured in 
a detector far away from the surface. By varying the delay between the pump and probe 
pulses, the intensity of the 2PPE signal as a function of the delay will reflect the evolution 
of the population of the state n. From this information it is then possible to extract the 
lifetime of the state n. 

Lifetimes of image states for different metal surfaces have been calculated in the frame- 
work of the self-energy formalism in the GW approximation 0, 0]. The results obtained 
happen to be in rather good agreement with experiments. For example, the calculated 
lifetime of the first three image potential states of the Cu(100) surface are rJ heor = 30 fs, 



r theor = 132 fs, and r| heor = 367fsj^|, while the experimental values are rf xp = 40 ± 6fs, 
r% xp = HOilOfs, and r^ xp = 300 ± 15 fs 61 , while for the first image state of Cu(lll]surface 
the theoretical lifetime is r^ heor = 17.5 fs^f and the experimental is r^ xp = 15 ± 5fs ^]. 

In the present paper we present calculations on energy-resolved 2PPE spectra of Cu(100) 
and Cu(lll) surfaces. Our theoretical approach is based on the propagation of an electron 
wave-packet in an one-dimensional model potential^! simulating the Copper surface. The 
finite electronic lifetimes of the image potential states are taken into account by an empirical 
self-energy term. The photoemission spectra are obtained through a geometrical scheme 
described in Ref. [ll|. Note that in this approach the potential is fixed during the whole 
simulations, i.e., the change of the electronic screening during the excitation is not taken 
into account. This could be important for short laser excitations as done in TR2PPE. 

To assert the relevance of this approximation we also performed time-dependent simula- 
tions of the Copper surface, but including the change of the Hartree potential due to the 
electronic excitations. 



METHOD 



To model the Copper surface we used a one-dimensional slab of 48 Copper layers, sur- 
rounded by 290 a.u. of vacuum on each side. The large portion of vacuum is necessary in 
order to describe the four first image potential wave functions (see Fig. [Q) . The bulk and 
image potentials are modeled by the ID potential model of Ref. Q]. This model reproduces 
the position and the width of the energy gap, as well as the energy of the surface state and of 
the first image potential state. Furthermore, it provides a good description of the electronic 

. In this model the number of bulk states is 
discrete for each value of q\\, and is related to the number of layers in the simulation. We do 
not believe that this is a major limitation, essentially because all the relevant physics occurs 
close to E-p, while the remaining states mainly act as a polarizable background. The elec- 
trons are allowed to move and interact only in the z-direction, perpendicular to the surface, 
and we assume a parabolic dispersion in the parallel direction. Within this approximation, 
the Hamiltonian describing such system can be written as 

Hq — ~~1~2 + "^model • (3) 
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This operator is discretized in real-space using a grid-spacing of 0.2 a.u, which is sufficiently 
small to allow for a proper description of the relevant electronic states. 

At t = we assume the system to be in the ground-state of the Hamiltonian (j^J. We 
then propagate this state with the new Hamiltonian 

H(t) = H + iwfo t) + t(z) , (4) 

where ■0i aS er(^, t) describes a laser polarized in the z direction that reads, in the length gauge, 
v\aser(z,t) = zE(t). The electric field is composed of a pump and a probe pulses, 

E(t) = E pump (t) cos(u pump t) + E prohc (t) cos(u prohc t) . (5) 

the functions -E pU m P / P robe(£) are envelope functions of the form 



Ei{t) = Ei cos 2 
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(t-ti 



Q{t-ti + <Til2)Q(-t + U + <Til2), (6) 



where ex, is the width of the pulse and tj the center of the pulse. For the pump pulse, t pump 
is simply o" pump /2, while for the probe pulse t pro be = + Cpump./S, with £ D the delay time 
between the pump and probe lasers. In order to take into account the finite lifetime of 
the image potential states we add to the Hamiltonian (JHJ) a time-independent self-energy 
operator of the form: 

N 

E(z) = -i J2 T k\n k ) (n k \ (7) 
fc=i 

where \rii) are the image potential states and are the inverse experimental lifetimes of 
each state taken from Ref. Note that H(0) = H . The time-propagation is performed in 
real-time following the method of Ref.jllj]. This propagation scheme preserves unitarity for a 
hermitian Hamiltonian, and has been proven very robust and stable in diverse applications. 



ENERGY-RESOLVED SPECTRA 



In order to obtain the energy-resolved photo-electron spectra we follow the technique 
of Ref. [lfll| and divide the simulation box in two regions: one region, A, containing the 
slab and the region where the first image potential states lie, and another, B, defined as 
the complement of A (see Fig. EJ). In region B the electrons are considered free outgoing 
particles, and are treated in momentum space. The separation between the two regions is 
achieved through a smooth masking function, M, defined as one in the interior of region A 
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FIG. 2: Typical shape of the mask function, M, used for dividing the space in two regions (see 
text for details). For illustrative purposes we also plot the densities of the first three image states 
of Cu(100). 



and zero outside. The method consists in evolving the wave-packet in the first region with the 
Hamiltonian (|%|). and mask the orbitals at each time-step. The electrons that "leave" region 
A are then treated as free-particles and accumulated in momentum space. In mathematical 
terms, the evolution is performed in the following way: We start by propagating the wave- 
functions in regions A and B using 

V>feCM + At) = M(z)exp(-iHAt)'il£(z,t) (8) 
i$(p,t + At) = exp | -i [P - ^ (t)] At j ^(p, t) + ffi(p, t + At), 

where i/)£(p,t + At) is the Fourier transform of the part of the wave function if)£ that left 
region A during the time step At, i.e, 

j£(p,t + At) = Jdz exp(ipz)[l - M(z)} exp(-iHAt)^(z,t) , (9) 

where p denotes the momentum. If we wait long enough after that the laser has been turned 
off, region B will contain those electrons that were ionized. The photo-electron spectra is 
then identified with: 

JV 



P(V2^E) = £ \^(p,t - oo)| 2 (10) 
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FIG. 3: Photoelectron energy resolved spectra of a Cu(100) (left panel) and a Cu(lll) (right panel) 
surfaces for both zero and finite pump-probe delay. 



where iV is the total number of electrons. This method can be derived from the interpretation 
of the Wigner transform of the one-body density matrix as a probability density [10]. 

In Fig. El we depict the energy resolved spectra of a Cu(100) and a Cu(lll) surfaces for 
both zero and finite pump-probe delay. For the (100) surface, the binding energies of the 
first image potential states are -0.57, -0.18, and -0.08eV (see Fig. |TJ), while for the (111) 
surface the first states are at -0.82, -0.22 and -0.009 eV. Note that the second and third image 
states of the Cu(lll) surface are resonances. Furthermore, the laser parameters used in these 
calculations were, for the (100) surface, hu pnmp = 4.7 eV, %u) pTO he = 1.8 eV, o" pump = 95 fs, and 
Cprobe = 54 fs, and for the (111) surface, hoj pumv = 5.1 eV, hu pTO h e = 1.8 eV, cr pump = 87 fs, 
and (T pro be = 50 fs. From these results we observe that the peaks corresponding to the first 
image potential state are obtained at the correct energy position. The signal due to the 
image states of higher quantum number cannot be distinguished due to the small difference 
in energy that cannot be resolved with our grid resolution. We also note that, for zero delay, 
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FIG. 4: Left: Population of the image states n = 1,2, and 3 as function of time. The values 
of rj = i 2 3 appearing in Eq. ((7J) were taken from Ref. |3j, and the maximum of each curve was 
normalized to one. Right: 2PPE intensity for the n = 1,2, and 3 image states of Cu(100) as 
function of the pump-probe delay. Reproduced from Ref. [g] 



the broadening of the peak relative to the first image potential state is larger in the case of 
the Cu(100) surface than in the Cu(lll). This is due to the fact that in the former case, the 
probe pulse can 'pump' low lying energy electrons from the bulk into energy levels near the 
gap which can be then ejected from the surface by the pump pulse. As the Cu(lll) surface 
has a different electronic structure, such situation does not occur[3|. Note that Cu(100) 
shows an intrinsic surface state resonance, while Cu(lll) has a surface state below the Fermi 
level. In the last case, the dominant transition happens between the occupied surface state 
and the first image potential state. This can be rationalized in terms of the larger spacial 
overlap between surface and image states as compared to bulk states. 



TIME-RESOLVED SPECTRA 



The typical energy width of the laser pulses used in TR-2PPE varies between about 10 
and 30 meV. With this energy resolution it is possible to excite separately each one of the 
first three image potential states of Cu(100). By following the TR-2PPE intensity in function 
of the pump-probe delay it is then possible to follow the evolution of the population of these 
states. In Fig0]we present calculations of the populations of the first three image states of 
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FIG. 5: Population of the first image potential state of Cu(100) calculated with and without the 
term of Eq. ljlljl for a 29 fs laser pulse of energy Tito = 4.2 eV. In both cases the maximum value of 
the population has been normalized to one. The difference in the lifetime for the two cases is 2% 



Cu(100). The curves were obtained by propagating the system with the Hamiltonian (J1J). 
For comparison we also show the experimental data from Ref. jf|. Clearly, there is a good 
qualitative agreement between the theoretical and experimental curves. 

In these calculations, the electron-electron interaction was taken into account by the 
model potential of Ref. |2j. However, the variation of this term with time was mostly 
neglected. Part of the time- dependence is taken indirectly into account through the self- 
energy operator, but one can question how reasonable this approximation really is. In order 
to answer to this concern, we performed simulations in which we allowed the classical part 
of the electron-electron interaction (i.e. the Hartree term) to change with time. This was 
achieved by adding to the Hamiltonian (jlj) a new term, 5i>Hartree(<z, t), that corresponds to 
the variation of the Hartree potential due to the change of electronic density in the direction 
perpendicular to the surface, 

— = -4n6p(z, t) . (11) 

where 5p(z, t) = p(z, t) —p(z, 0). Exchange and correlation effects were treated as before, i.e. 
through the self-energy term (JZj) used to simulate the lifetime of the states. The evolution 
of the population of the first image state with and without this term is shown in Fig 03 We 
observe only a minor shift of the peak to lower times from 24fs to 21.3fs and the lifetime is 



reduced by around 2% when the term of Eq. (fTTj) is included. This result validates previous 
studies that have neglected this effect Q]. 



CONCLUSIONS 

We have performed dynamical simulations in order to calculate energy-resolved photoe- 
mission spectra of the Cu(100) and Cu(lll) metal surfaces. Using the same technique we 
also obtained time-resolved photoemission spectra for individual excitations of the first three 
image potential states of Cu(100). This technique is quite general, and can be used to gain 
insight on the dynamics of image states with large quantum number and when more than one 
eigenstate is excited coherently (the situation in quantum beat spectroscopy). Furthermore, 
we studied the influence of the variation of the Hartree potential during the excitation. This 
effects turns out to be small (around 2%). 

Acknowledgments 

This work was supported by the Basque Country University, DGES and European Com- 
munity under the research training network NANOPHASE (HPRN-CT-2000-00167). We 
wish to thank H. Appel, E. K. U. Gross and P. M. Echenique, for enlightening discussions 
and to the high-performance computing facilities of CEPBA where most of the calculations 
were done. 



[1] P. M. Echenique and J. B. Pendry, J. Phys. C 11, 2065 (1978); P. M. Echenique and J. B. 

Pendry, Prog. Surf. Sci. 32, 111 (1990); P. M. Echenique, J. M. Pitarke, E. V. Chulkov, V. M. 

Silkin J. El. Spec. Rel. Phen. 126, 163 (2002). 
[2] E. V. Chulkov, V. M. Silkin and P. M. Echenique, Surf. Sci. 437, 330 (1999). 
[3] P. M. Echenique, J. M. Pitarke, E. V. Chulkov and A. Rubio, Chemical Physics 251, 1 (2000); 
[4] Chemical Physics 251, 1 (2000); T. Fauster, W. Steinmann, in: P. Halevi (Ed.), Photonic 

Probes of Surfaces, Electromagnetic Waves: Recent Developments in Research Vol.2, Elsevier, 

Amsterdam, 1995. 
[5] K. Giesen et al, Phys. Rev. Lett. 55, 300 (1985). 

10 



[6] H. Heifer, I. L. Shumay, C. ReuB, U. Thomann, W. Wallauer and T. Fauster, Science 277, 
1480 (1997). 

[7] T. Klamroth, P. Saalfrank and U. H6fer, Phys. Rev. B 64, 035420 (2001). 

[8] J. Osma, I. Sarria, E. V. Chulkov, J. M. Pitarke and P. M. Echenique Phys. Rev. B 59, 10591 
(1999); E. V. Chulkov, J. Osma, I. Sarria, V. M. Silkin and J. M. Pitarke Surf. Sci. 433, 882 
(1999); I. Sarria, J. Osma, E. V. Chulkov, J. M. Pitarke and P. M. Echenique Phys. Rev. B 
60, 11795 (1999); E. V. Chulkov, I. Sarria, V. M. Silkin, J. M. Pitarke and P. M. Echenique 
Phys. Rev. Lett. 80, 4947 (1998) 

[9] J. D. McNeil, N. H. Ge, C. M. Wong, R. E. Jordan, C. B. Harris, Phys. Rev. Lett. 79, 4645 
(1997). 

[10] D. Varsano, M. A. L. Marques, H. Appel, E. K. U. Gross and A. Rubio, to be published. 
[11] M. A. L. Marques, A. Castro, G. F. Bertsch and A. Rubio, Comput. Phys. Commun. in press 
(2002). 

[12] For a schematical projection of the bulk band structure onto the (111) and (100) surfaces of 

Cu see, for instance, Ref. 3] pag. 21. 
[13] This expression corresponds to assuming an infinite repulsive barrier at the surface. In the 

more general case we have V(z) = — u z e _ z . where z- lxa is the position of the image plane 

(that depends on the electronic and structural properties of the surface), and e is the static 

dielectric constant of the surface. 



11 



